VARIATIONAL APPROACH TO COMPLICATED SIMILARITY 
SOLUTIONS OF HIGHER-ORDER NONLINEAR PDES. II 



V.A. GALAKTIONOV, E. MITIDIERI, AND S.I. POHOZAEV 

Abstract. This paper continues the study began in [TTJ [T^] of the Cauchy problem 
for {x, t) e X for three higher-order degenerate quasihnear partial differential 
equations (PDEs), as basic models, 

ut = (-I)"+iA'"(|m|"m) + 
utt = (-I)™+iA™(|u|"u) + 

where n > is a fixed exponent and A™ is the (to > 2)th iteration of the Laplacian. 
A diverse class of degenerate PDEs from various areas of applications of three types: 
parabolic, hyperbolic, and nonlinear dispersion, is dealt with. General local, global, and 
blow-up features of such PDEs are studied on the basis of their blow-up similarity or 
travelling wave (for the last one) solutions. 

In [HI [12], Lusternik-Schnirel'man category theory of variational calculus and fiber- 
ing methods were applied. The case m = 2 and n > was studied in greater detail 
analytically and numerically. Here, more attention is paid to a combination of a Carte- 
sian approximation and fibering to get new compactly supported similarity patterns. 
Using numerics, such compactly supported solutions constructed for m = 3 and for 
higher orders. The "smother" case of negative n < is included, with a typical "fast 
diffusion-absorption" parabolic PDE: 

Ut = (-l)™+iA™(|u|"u) - Iwr-a, where n G (-1,0), 

which admits finite-time extinction rather than blow-up. Finally, a homotopy approach is 
developed for some kind of classification of various patterns obtained by variational and 
other methods. Using a variety of analytic, variational, qualitative, and numerical meth- 
ods allows to justify that the above PDEs admit an infinite countable set of countable 
families of compactly supported blow-up (extinction) or travelling wave solutions. 



1. Introduction: higher-order blow-up and compacton models 

A general physical and PDE motivation of the present research can be found in flT\ [T2] , 
together with basic history and related key references, so we just briefly comment on where 
quasilinear elliptic problems under consideration are coming from. 



Date: January 14, 2013. 

1991 Mathematics Subject Classification. 35K55, 35K40, 35K65. 

Key words and phrases. Quasilinear DDEs, non-Lipschitz terms, similarity solutions, blow-up, extinc- 
tion, compactons, Lusternik-Schnirel'man category, fibering method. 

1 



1.1. (I) Combustion- type models with blow-up. Our first model is a quasilinear 
degenerate 2mth-order parabolic equation of the reaction-diffusion (combustion) type: 

(1.1) ut = (-1)''"+1A'"(|m|"m) + l^rw in X R+, 

where n > is a fixed exponent, m > 2 is integer, and A denotes the Laplace operator 
in M^. Physical, mathematical, and blow-up history of (11. ip for the standard classic case 
171 = 1 and m > 2 is explained in [1111121 § 1]. Consider regional blow-up solutions of (II. ip 

(1.2) us{x,t) = {T -t)-^f{x) in M^x(0,T) 

in separable variables, where T > is the blow-up time. Then the similarity blow-up 
profile f = f{x) solves a quasilinear elliptic equation of the form 

(1.3) (_i)-+iA-(|/|"/) + |/|"/ = i/ in R^. 

This reduces to the following semilinear equation with a non-Lipschitz nonlinearity: 





F=\frf =^ 
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in the nonlinear term yields 
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According to (11.21) . the elliptic problems (11. 4p and the ODE (II. 5p for = 1 are responsible 
for the possible "geometrical shapes" of regional blow-up described by the higher-order 
combustion model (II. ip . 

Remark: relation to ODEs from extended KPP theory. There exists vast mathe- 
matical literature, starting essentially from the 1980s, devoted to the fourth-order ODEs 
(looking rather analogously to that in (II. 5p for m = 2) 

(1.6) F^^) = + F -F^ in M, 

where /3 > is a parameter. This ODE also admits a complicated set of solutions with 
various classes of patterns and even with chaotic features. We refer to Peletier-Troy's 
book [26] for the most diverse account, as well as to papers [TH [31], where a detailed 
and advanced solution description for (II. 6p is obtained by combination of variational 
and homotopy theory. Regardless their rather similar forms, the ODEs (II. 5p belong 
to a completely different class of equations with non-coercive operators, unlike in (II. 6p . 
Therefore, direct homotopy approaches and several others, that used to be rather effective 
for (II. 6p . fail in principle for (II. 5p . In this sense, (II. 5p is similar to the cubic ODE to be 
studied in Section |6l 

(1.7) 



m 



of course, excluding complicated oscillatory behaviour at finite interfaces, which are obvi- 
ously nonexistent for analytic nonlinearities. However, we claim that the sets of solutions 
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of (11 -Sp . with m = 2, and of (II .Tp are equivalent, though, not having a rigorous proof, 
we will devote some efforts to a homotopy approach connecting solutions of such smooth 
(analytic) and non-smooth ODEs. Thus, though going to develop homotopy approaches 
for classifying solutions of (ll.Sp (Section H]), our main tool to describe countable families 
of solutions {Fi} is a combination of Lusternik-Schnirel'man category-genus theory [21] 
and the fibering method [271 [2H] ■ 

Thus, we show that ODEs (II. 5p . as well as the PDE (II ■4p . admit infinitely many count- 
able families of compactly supported solutions, and the whole solution set exhibits certain 
chaotic properties. Our analysis will be based on a combination of analytic (variational 
and others), numerical, and some more formal techniques. Explaining existence, multiplic- 
ity, and asymptotics for the nonlinear problems involved, we state and leave several open 
difficult mathematical problems. Meantime, let us characterize other models involved. 

1.2. (II) Regional blow-up in quasilinear hyperbolic equations. Consider next 
the 2mth-order hyperbolic counterpart of (II. ip . 

(1.8) Utt = (-1)'"+^A"(|m|"m) + in x M+. 

The blow-up solutions take a similar form with a different exponent — - instead of — -: 

(1.9) 7/s(x,t) = (T-t)-^/(a;), 

where F = | /!"■/, after scaling, solves the same elliptic equation (II. 4p . 

1.3. (Ill) Nonlinear dispersion equations and compactons. Such rather unusual 
PDEs in A^-dimensions (the origin is integrable PDEs and other areas) take the form 

(1.10) ni = (-l)'"+i[A™(|nrn)],, + (|Mrn),, in x M+, 

where the right-hand side is the derivative D^^ of that in the parabolic counterpart (II. ip . 
Then the elliptic problem (II. 4p occurs when studying travelling wave (TW) solutions of 
(ll.lOp . Note that, as being PDEs with nonlinear dispersion mechanism, (ll.lOp and other 
NDEs listed below admit shock waves and other discontinuous solutions. Here, we study 
smooth solutions and do not touch difficult entropy-hke approaches for such shock and 
rarefaction waves and refer to [T^ § 4.2] and [13] for an account to such phenomena. 

Thus, for the PDE (ll.lOp . looking for a TW compacton (i.e., a solution having all the 
time compact support; see key references in pj, § 1]) moving in the xi-direction only, 

(1.11) Uc{x,t) = f{yi,X2, ...,xn), yi=xi-Xt, where A = -i, 

we obtain on integration in yi the elliptic problem (II. 4p . Analogously, for the higher-order 
evolution extension of nonlinear dispersion PDEs, 

D^u = D^^ [(-1)'"+iA'"(|u|"m) + in x M+ {k > 2), 

to get the same PDE (II. 4p . the compacton (II. lip demands the following wave speed: 



1.4. (IV) PDEs with "fast diffusion" operators: parabolic, Schrodinger, hy- 
perbolic, and nonlinear dispersion models. This is about negative exponents: 

(1.12) nG(-l,0), 

which generate other types of elliptic equations of interest. To connect such problems 
with typical models of diffusion-absorption type, consider the following parabolic PDE: 

(1.13) ut = (-1)'"+^A'^(|m|"m) - in x M+ {-1< n < 0), 

with the strong non-Lipschitz at m = absorption term. It is well known that such 
PDEs describe finite-time extinction phenomenon instead of blow-up. See [161 Ch. 4,5] 
for m = 1 and [U [32] for m > 2 for necessary references and history of strong absorption 
phenomena. Therefore, the similarity solution takes an analogous to fll.2p form with the 
positive exponent — - > 0, so that u{x,T~) = 0, while / solves a similar elliptic equation 

(cf. m) 

(1.14) (-l)™+iA™(|/r/) - l/r/ = i / in M^. 

By the scaling as in (11. 4p (recall that n < 0), we eventually obtain the semilinear elliptic 
problem with a sufficiently smooth nonlinearity: 



;i.i5) 



-l)™+iA™F - F + F "F = in M^, where a = — ^ > 0. 

' ' n+l 



The nonlinearity is now at F = 0, so the solutions are classic. For instance, for m = 2 
and n = — |, we obtain equation with a cubic analytic nonlinearity: 

(1.16) A'^F = -F + F^ in R^. 

For = 1, this is the ODE (11. 7p . Indeed, these equations do not admit solutions with 
finite interfaces and exhibit exponentially decaying oscillatory behaviour at infinity. We 
show that the total set of such "effectively" spatially localized patterns well matches those 
ones for ?7, > always having finite interfaces. 

As a connection to another classic PDE area and applications, let us note that, for 
m = 1, we obtain the classic second-order case of the ground state equation [6J 

(1.17) AF-F + F^ = in M^. 

This elliptic problem is key in blow-up analysis of the critical nonlinear Schrodinger 
equation (NLSE) 

iut = —Au — lul'^^^u, p = po = l + j^, N = 2, and ^(x, t) = e^*F(x); 

see Merle-Raphael [22]-[2l] as a guide. The solution F of (11.170 is strictly positive (with 
exponential decay at infinity) and is unique up to translations, while the ground state Fq 
for (I1.16P is oscillatory at infinity, to say nothing about a huge variety of other, Lusternik- 
Schnirel'man or not, solutions. Thus, (I1.16P is the ground state equation for the fourth- 
order NLSE 

iut = A'^u-\u\P-^u in R^ x R+, where p = pQ = 1 + {N = 4). 



Analogous similarity analysis is performed for the corresponding "fast diffusion" hyper- 
bolic equation (the extinction patterns are given by fll.9p ) 

(1.18) utt = (-l)™+'A"(|MrM) - in x M+ (-K n < 0), 
and for the nonlinear "fast dispersion" PDE 

(1.19) Mi = (-l)™+i[A"(|MrM)],, -(1^1"^),, in M^xM+ (-K n < 0), 
where the moving TW profiles are as in (II. lip . 

1.5. Main goals and connections with previous results. It turns out that such 
profiles F solving (11. 4p have rather complicated local and global structure. The study of 
equations (ll.4p and (II. 5p was began in [11], where the following goals were posed: 

(i) Problem "Blow-up": proving finite-time blow-up in the parabolic (and hyperbolic) 
PDEs under consideration [H] § 2]; 

(ii) Problem "Existence and Multiplicity": existence and multiplicity for elliptic PDEs 
([11]) and the ODEs ([L5]) for m > 2 § 3]; 

(ill) Problem "Oscillations" : the generic structure of oscillatory solutions of (II. 5p near 
interfaces for arbitrary m > 2 [111 § 4]; and 

(iv) Problem "Numerics" : numerical study of all families of F{x) for m = 2 [TT] § 5]. 

The non-Lipschitz problem (ll.ip possesses so complicated set of admissible compactly 
supported solutions (note that Lusternik-Schnirel'man category theory detects only a 
single countable subset), that using effective MatLab (or other similar or advanced) nu- 
merical techniques for classifying the critical points becomes an unavoidable tool of any 
analytic-numerical approach, which cannot be dispensed with at all. We recall that in [TTl 
§ 3] the identification of Lusternik-Schnirel'man sequence of critical values was confirmed 
numerically only (and for m = 2 essentially). 

Therefore, in the present paper, in Section [3], we begin by continuing achieving the goal 

(iv) for the sixth-order case m = 3: 

(iv') Problem "Numerics" : numerical study of all families of F[x) for m > 3 (Section |3]). 

In addition, we also aim new targets: 
(ii') Problem "Existence and Multiplicity": using variational Lusternik-Schnirel'man ap- 
proach and fibering with an auxiliary Cartesian approximation of critical points (Section 
E]); 

(v) Problem "Fast diffusion": n G (—1,0), where smoother elliptic problems (ll.lSp occur 
(extinction in Sections |5] and existence- multiplicity in Section [6]). 

Finally, for both the non-smooth (11. 4p and smooth ( ll.lSp problems, we pose: 

(vi) Problem "Sturm Index" : classification of various patterns according to their spatial 
shape for both (ll.4p and (ll.lSp (Sections H] [6l and [7]). [For m = 1, this is governed by 
classic Sturm's First Theorem on zero sets.] 

Thus, we are introducing three classes, (I) , (II) , (III) , of nonlinear higher-order PDEs 
in X M , including similar fast diffusion ones (IV). These are representatives of PDEs 
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of three different types. However, it will be shown that these exhibit quite analogous 
evolution features (if necessary, up to replacing blow-up by moving travelling waves or 
extinction behaviour), and coinciding complicated countable sets of evolution patterns. 
This reveals an exiting feature of a certain unified principle of singularity formation phe- 
nomena in general nonlinear PDE theory, which we would like to believe in, but which is 
very difficult to justify more rigorously. 

1.6. On extensions to essentially quasilinear equations. The three-fold unity of 
PDE classes (I)— (III) is available for other types of nonlinearities. In Section |8l we 
briefly discuss the following classes of equations with fourth-order p-Laplacian operators 
(here n = p — 2 > 0): 

(I) Ut = — A(| A-uI^Am) + \u\"'u (parabolic), 

(1.20) (II) utt = -A(|Am|"Am) + (hyperbolic), 
(III) Ut = -[A{\Au\"Au)]^^ + {\u\^u),, (NDE). 

It turns out that these equations admit similar blow-up or compacton (for the NDE (III)) 
solutions that are governed by variational elliptic problems with similar countable variety 
of oscillatory compactly supported solutions. As a first step, an approach to blow-up of 
solutions of the parabolic equation (11.201) for = 1 and some other related results on 
similarity solutions can be found in [9]. 

Thus, these cases are more difficult and are essentially quasilinear, since the resulting 
elliptic problems cannot be reduced to semz/m ear equations such as (11. 5p . 

1.7. Towards non-variational problems: branching. Principally more difficult prob- 
lems occur under a slight change of the lower-order nonlinearity in fll.20p : 

(1.21) h-). \u\^u, with p > n. 

This leads to non-variational elliptic problems, to be also briefly discussed in Section [8] 
using the idea of branching of proper solutions aXp = n from the similarity proflles studied 
above. Proving existence of countable sets of solutions for such non-potential operators 
reveals a number of open problems of a higher level of complexity. 

2. Other families of patterns: Cartesian approximation and fibering 

Application of the Lusternik-Schnirel'man category theory to constructing a countable 
family of solutions of (11.51) is explained in [llj § 3]. This allowed us to detect the so-called 
basic family of patterns {-F]}, which has been shown for A^ = 1 in a number of flgures in 
[m § 4]. For convenience, we restate the Lusternik-Schnirel'man/fibering result in [TTl 
§ 3.2]: 



Proposition 2.1. The elliptic problem ^1.4\ has at least a countable set of different 



solutions denoted by {Fi, I > 0}, each one Fi obtained as a critical point of the functional 
(2.1) E{F) = -\I \D-F\^ + \I F'-\l \F\^ /3 = ^ G (1, 2), n > 0, 



in Wq'^IBpi) in a ball Br with a sufficiently large radius R = R{1) > 0. 

2.1. Basic computations. We next develop approaches for obtaining other patterns, 
which are not detected in Proposition 12 . 1 1 by Lusternik-Schnirel'man and classical fibering 
techniques. 

In order to construct other families of solutions (see Section [3] for illustrations of those 
for m = 3), we need an auxiliary approximation of patterns. Namely, we first perform 
the Cartesian decomposition 

(2.2) F = h + w, 
where h G (^ij) 

is a smooth "step-like function" that takes the equilibrium values ±1 
and on some disjoint subsets of Br (with a smooth connections in between). Sufficiently 
close to the boundary points, we always have h{y) = 0. For instance, in one dimension 
for getting the patterns in Figured! we take h{y) as a smooth approximation of the step 
function, which takes values ±1 and on the intervals of oscillations of the solution about 
these equilibria. 

In other words, we are going to perform the radial fibering not about the origin but 
about the non-trivial point h, which plays the role of an initial approximation of the pat- 
tern that we are interested in. Obviously, the choice of such /I's is of principal importance, 
which thus should be done very carefully. 

Substituting (12.21) into the functional yields the new one, 

E{w) = E{h + w) = -- f \D'^h\^ + - [ h^ + Lo{h)w 

(2.3) ^^'^ 

-Ho{w)-- / \h + w\^, 

P J Br 

where by Lq we denote the linear functional 



Lo{h)w = - j D'^h ■ D'^w + I hw. 

J Br J Br 



We next apply the fibering approach by setting, as usual, 
(2.4) w = r{v)v, V G Hq, whence 

H{r, v) = E{h + r{v)v) = -- [ \D'^h\'^ + - / + r Lqv 



2 J Br 2 



1.1/,,. 

In order to find the absolute minimum point, we need to solve the scalar equation H'^ = 0, 

(2.6) Hr{r,v)=r- \h + rvf''^{h + rv)v + Lo{h)v = 0. 

Jbr 



m=2, n=1 : pattern F 




m=2, n=1 : pattern F 



■6,2,+4.1,-2,1,+2 



-30 -20 -10 10 20 30 



(a) (7 = {+6, 2, +2, 2, +6} 



m=2, n=1 : pattern F 



■2,2,+4,2,+2,1,-4 




10 20 




(b) a = {+6, 2, +4, 1,-2, 1 + 2} 



m=2, n=1: pattern F^^ ^ _^ ^ 




(c) o = {+2, 2, +4, 2, +2, 1, -4} (d) a = {+6, 3, -4, 2, -6} 

Figure 1. Various patterns for m = 2 and n = 1. 



For /i = 0, this coincides with the standard equation derived in [HI § 3], and has three 
roots, ro(f ) = and 



(2.7) 



r± = r±[v) 



which are positive and negative respectively. For h ^ 0, these roots exist and are shghtly 
deformed for sufficiently small h. For large h, one of the roots r±{v) may disappear, 
and at this instance the resulting functional {q.v. below) may loose its smoothness. To 
distinguish the roots, we observe that 



(2i 



r„(t>) < always exists and smooth if J(/i, v) = — '^hv + Lo(/i)f > 0, 

J Br 

r+(t>) > always exists and smooth if J{h,v) = — \hf~'^hv + Lo{h)v < 0. 



Thus, calculating the extremum point r±{v) (when exists) from fl2.6p and substituting 
into fl2.5p yields the new functional 

(2.9) H{v) = H{r±{v),v) on Hq, being even, since r±{—v) = —r±{v). 

Therefore, if fl2.9p is smooth on an appropriate branch r±, this gives a set of critical point 
{vk}, as above. Moreover, in a neighbourhood of any critical points Vk satisfying (12. 8p . 
i.e., J{h,Vk) > or < 0, the corresponding branches r±{v) are smooth for v ^ Vk (and 
hence along a minimizing sequence {vl} — )■ Vk), so that (12. 9p is sufficiently regular. Even 
in the delicate case, when J{h,Vk) = and r = is an inflection point of Hr{r,Vk), i.e., 

Hrr{r,Vk) = 0, 

one can choose a smooth, existing, and "stable" branch r = ro(f) for v ^ Vk along a 
suitable minimizing sequence. 

This provides us with a finite number of critical points associated with the category of 
Hq. As we have seen, all these critical points of (12. 3p such as 

— sign {Lo(/i)w} = sign {Lo(/i)f } = ±1 (or 0) 

can be obtained on the branch r±{v) (or ro(t')). 

Actually, using this mixture (12. 2p of the Cartesian and spherical fibering decomposition 
of the functional space, we are interested, mainly, in the first critical point, which is defined 
via the absolute infimum of the functional (12. 9p (roughly speaking, in the case of genus 
1). With a choice of a sufficiently "large" approximating function h, this first pattern 
will be different from other basic patterns constructed above for /i = 0. Indeed, this first 
pattern is characterized by the condition of the "minimal deviation" from h{y), while, 
e.g., Fq corresponds to the minimal deviation from /i = 0, so that these cannot coincide 
if h is large enough and has a proper shape concentrating about equilibria ±1 and 0. 

Figure [2] illustrates such a statement and shows a typical Cartesian approximation h, 
which is necessary to detect the patterns F+4^_4^+2,_2,+2- Obviously, then the absolute 
extremum of H{v) cannot be attained at already known critical point Fq given by the 
dashed line, which is characterized by a much larger deviation from the fixed h that is 
given by a boldface line (it should be slightly smoothed at corner points). 

Therefore, the main result in [TT| § 3] , such as Proposition 12. ![ remains true for any 
sufficiently regular initial approximation h. Of course, some of the critical points Fi{y] h) 
with I ^ 1 may coincide with already known basic patterns Fi, but, in fact, we are 
interested in the first critical value and point, which thus give F that has the minimal 
deviation from h, and must be different from F/'s. Obviously, for approximations h that 
are far away from 0, the first pattern F{y] h) obtained by using the sets T-Lq of arbitrary 
category, p > 1, cannot coincide with the first basic patterns Fi{y), which are sufficiently 
small and have a specific and different geometric structure. 

The actual and most general rigorous "optimal" choice and characterization of such 
suitable approximations h (possibly a sequence of such {hk}) remains an open problem, 
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Hy) 


Hy) 


FJv) -^+4-4,+2-2,+2(Z/) 

--4 LX A- 















Figure 2. For getting F4.4^_4^+2,-2,+2) the Cartesian approximation h in 
should be chosen properly. 

though we have got a convincing experience in understanding of such patterns, in par- 
ticular, using numerical experiments and some analytic estimates; see related comments 
below. 

2.2. Some asymptotic analysis. It is easy to show that, asymptotically, for sufficiently 
"spatially wide" patterns, the Cartesian-spherical fibering (12. 2p . (12. 4p provides us with 
families of patterns that are different from basic ones {Fi}. 
For instance, in Figure [3] we compare the patterns 

F+2fc (the dashed line) and -F+2,2,+2,2,...,+2,2,+2 for large k = 10, i.e., 
F+20, where 



Cf = 2.9398... 



and 



where 



Cf = 2.7197... 



-P+2,2,+2,2,+2,2,+2,2,+2,2,+2,2,+2,2,+2, 

Here cf stand for the corresponding critical values of the functional obtained after fibering 

[mini §3] 



(2.10) 



Cf = H{v) 



\F\P 



Be 



,3/2 



(/5 



n+2 \ 
n+l)' 



Thus, these two patterns are clearly recognized by their different critical values cf 
indicated. For -F+20, the corresponding functions h{x) 1 on (—40,40) is shown by the 
boldface line. It is seen that the global minimum of the functional (12. 9 p for such h cannot 
he attained on any profile from the basic family {-F]}, because the total deviation becomes 
huge in comparison with the almost periodic deviation achieved via F+2o- In this case, 
the minimum is attained on the profile Fj^2k having a completely different geometry. 
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m=2, n=1 : profiles F^^o and F^2,2,.2,2,.2,2,.2,2,.2,2,.2,2,.2,2,.2,2 

1.6 
1.4 



1.2 

F(y)i 

0.8 
0.6 
0.4 
0.2 


-40 -30 -20 -10 10 20 30 40 

y 

Figure 3. Two patterns, (dashed line) obtained by Cartesian-spherical 

fibering and F+2,2,+2,2,+2,2,+2,2,+2,2,+2,2,+2,2,+2 corresponding to h = 0; m = 2, 
n = 1. 

For A; ^ 1, these observations can be fixed in a standard asymptotically rigorous 
manner, which we are not going to do here. 

2.3. The origin of countable sequences of solutions: a formal double fibering. 

Taking into account both changes (12.21) and (12 ■4p . we arrive at the functional 

(2.11) E{h,r,v) = E{h + r{v)v). 
The relative critical point of (I2.1ip are given by the system 

(E',{h,r,v) = 0, 

(2.12) lE'^{h,r,v) = 0, 

(E',{h,r,v) = 0. 

Of course, the first equation is just equivalent to the original one, since 

E'h = E', 

so that (I2.12P is a formal system comprising the spherical fibering in the {r, f }-variables, 
and the original equation. Let us see what kind of conclusions can be derived from this. 

The second equation is scalar and gives us necessary smooth branches (under certain 
hypotheses as above) 

(2.13) r = r*(/i, f) {r^{h, —v) = —r^{h,v)). 
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Then we arrive at the system 

, . fE'^{h + n{h,v)v) = 0, 

^ ^ \E'^{h + n{h,v)v) = 0. 

The first equation is difficult to handle (possibly more difficult than the original one). 
However, assume that this can be solved. In view of f l2.13p . gives an even dependence, 

(2.15) h = K{v) (Ki-v) = K{v)). 
Finally, we arrive at the even weakly continuous functional for v, 

(2.16) E,{v) = E{K{v)+r{K{v),v)v) in -Ho- 

This has a countabl^ set of critical point {vk}, which by fibering method [28], generate 
critical point of the original functional (12. lip . 

Therefore, eventually, as i? — oo, we obtain a countable set of necessary critical points 

{hk,Wk) = {{K{vk),r^{vk)vk), k = 0, 1,2, ... . 

The corresponding F-patterns are denoted by 

{Fi}. 

The actual general structure of such special solutions remains unclear and needs extra 
analysis. Currently, we know a little about this and present a few comments only. Using 
an analogy with the basic Lusternik-Schnirel'man patterns {Ei} obtained in fHX § 3] for 
h = 0, it may be expected that each E is composed from I > 1 copies of the "elementary" 
profile F+4, i.e., 

F,~{F+4,-F+4,...,(-l)'+'F+4}, 
with the obvious choice of the corresponding Cartesian approximation hi{y) that is con- 
centrated about equilibria at ±1 and in between. It is more likely that {Ei} includes 
other profiles of the {±F+4, ±F+4}-gluing (see [HI § 4] for definitions), or, in particular, 
can be composed from completely "non-oscillatory" profiles, i.e., 

Fi = E^2i- 

3. Problem "Numerics" : PATTERNS IN HIGHER-ORDER CASES, m > 3 

The main features of the pattern classification by their structure and computed critical 
values for m = 2 [11] can be extended to arbitrary m > 3 in the ODEs (11. 5p . so we 
perform this in less detail. 

In Figure m for the purpose of comparison, we show the first basic pattern Eo{y) for 
n = 1 in four main cases m = 1 (the only non-negative profile by the Maximum Principle 
known from the 1970s [31], [301 Ch. 4]), and m = 2, 3, 4. Next Figure |5] explains oscillatory 
properties of such Eo{y) close to the interface points. It turns out that, for m = 4, the 
solutions are most oscillatory, so it is convenient to use this case for further illustrations. 

In the log-scale, the zero structure of Eo{y) near interfaces is shown in Figure [6] for 
m = 2,3, and 4 {n = 1). For m = 4 and m = 3, this makes it possible to observe a 

^As usual, we mean compactly supported solutions in R^, i.e., R = +oo. 
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First pattern ^^=^^2 '^'"^'^ 





m=1 




m=2 




F(y)o 



Figure 4. The first solution Fo(y) of ([LS]), n = 1, for m = 1,2,3,4. 



First pattern F for m=1 ,2,3,4: oscillations near interface 




(a) scale 10 



First pattern F for m=1 ,2,3,4: oscillations enlarged 




(b) scale 10" 



FlGURE 5. Enlarged zero structure of the profile -Fo(y) for n = 1 from Figure 
HJ the linear scale. 



dozen of "nonlinear" oscillations that well correspond to the already known oscillatory 
component structure close to interfaces; see [llj § 4]. For the less oscillatory case m = 2, 
we observe 4 reliable oscillations up to 10^^°, which is our best accuracy achieved. 

The basic countable family satisfying approximate Sturm's property has the same topol- 
ogy as for m = 2 [TT, § 5], and we do not present such numerical illustrations. 
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m=2,3,4, n=1 : oscillations in log-scale 




y 

Figure 6. Behaviour of FQ(y) for n = 1, for m = 2,3,4; the log-scale. 



m=3, n=1 : patterns F^^^, k=1 ,2,3,4,5,6,7 
1.5 1 > > 1 > 




-30 -20 -10 10 20 30 



y 

Figure 7. The first seven patterns from the family {F^2k}', m = 3 and n = 1. 

In Figure [7] for m = 3 and n = 1, we show the first profiles from the family {F+2fc}, 
while Figure E] explains typical structures of F+2,fc,+2} for m = 4, n = 1. In Figure |9] for 
m = 4 and n = 1, we show the first profiles from the family {F+2,fe,-2}- 

14 



m=4, n=1 : seven symmetric patterns F^^ ^ for k=2,4,6 

1.5 



1 

F(y) 

0.5 







-20 -15 -10 -5 5 10 15 20 

y 

Figure 8. The first patterns from the family {F+2,fe,+2} of the {Fo,Fo}- 
interactions; m = 4 and n = 1. 

Finally, in Figure [IDl for comparison, we present a complicated pattern for m = 3 and 
4 (the bold line), n = 1, with the index 

(3.1) a = {-8, 3, +4, fc, -10, 1, +8, /, -12}. 

Both numerical experiments were performed starting with the same initial data. As a 
result, we obtain quite similar patterns, with the only difference that, in (13.11) . k = 1, 
I = 3 for m = 3, and for more oscillatory case m = 4, the number of zeros increase, so 
now A; = 3 and I = 5. 

4. Problem "Sturm Index" : AN HOMOTOPY CLASSIFICATION OF PATTERNS VIA 
e-REGULARIZATION (ANALYTIC-NUMERICAL APPROACH) 

4.1. Sturm index for second-order ODEs: in need to extension for m > 2. As 

we have mentioned, it is well-known that in the second-order case m = 1, solutions of 
ODE problems, even in the non-Lipschitz case, 

(4.1) F" = -F +\f\ in M 

obey Sturm's Theorem on zeros, which is a corollary of the Maximum Principle. Namely, 
concerning the problem (14. ip . each function Fi{y) from the basic family has precisely / 
isolated zeros (sign changes) and / + 1 non-degenerate extremum points. Therefore, the 
Sturmian index Ii = I of Fi, as the number of its "transversal zeros", uniquely specifies 

15 




m=4, n=1 : six "dipole" patterns F^^^ i _2k' .3,5,7 

1.5 I 1 . . 1 r-^ . 




y 

(b) zero structure on (—8, 8) 

Figure 9. The first patterns from tlie family {F+2,k,-2} of the {-Fq^Fq}- 
interactions, for m = 4 and n = 1: profiles (a), and zero structure around y = 



any of the basic patterns. This is also equivalent to the Morse index of the correspond- 
ing linearized operator. Moreover, the Lusternik-Schnirel'man minimax construction of 
critical points reveals this zero structure of the minimizers [211 P- 385] that is directly 
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m=3 and 4, n=1 : complex pattern F_j 




y 

Figure 10. A complicated pattern Ffj{y) for m = 3,4 and n = 1. 

associated with the category [3^, § 6.6], or genus [2T1 § 57], of the sets involved in the 
variational construction. 

For any m > 2, the Maximum Principle fails, and such a rigorous geometric classi- 
fication of basic patterns is no longer available in view of existence of oscillatory tails 
close to both interfaces. Roughly speaking, each profile obtained by the Lusternik- 
Schnirel'man/fibering approach has infinitely many zeros and extremum points, which 
makes it impossible to use the above simple geometric characteristics for classification of 
the patterns as for m = 1. 

Nevertheless, we claim that a Sturmian-type characterization of some (basic) patterns 
is possible for oscillatory solutions of higher-order ODEs. We will reveal how to attach the 
Sturmian index to solutions from the basic family {Fi] in the higher-order case, and also 
to other families. We consider the ODE case (11. 5p for = 1, though a similar approach 
applies to the radial elliptic setting in (II. 4p . as well as non-radial, where though it is not 
that well-presented and clear. 

We begin with description of higher-order equations admitting a rigorous Sturmian 
classification of patterns. 

4.2. 2mth-order equations with Sturmian ordering. Consider the following ODE 
with Dirichlet boundary conditions: 

(4.2) (_i)™+iF(2'n) + = in {-R, R), F = F' = ... = F^"^"^) = at y = ±R. 

This problem is also variational and admits a countable set of solutions {Fi}. Moreover, 
since the differential operator (— 1)™D^"^ is an iteration of the positive operator —D^ with 
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F*'''=F^, Dirichlet: first eight solutions, Tol=1Cr® 




Figure 11. The first eight patterns Fi{y) satisfying (fOj) for m = 2. 

the Maximum Principle, according to Elias [7] (see also applications to some nonlinear 
higher-order eigenvalue problems in [29l [2]). the following result holds: 



Proposition 4.1. The Ith-solution Fi{y) of the problem i \4-S^ for any I = 0, 1,2, has 
precisely I zeros and I + 1 extremum points on {—R, R). 



In Figures [H] and [121 we present the first solutions of (14. 2p for m = 2 and m = 4. 
According to Lusternik-Schnirel'man theory [2T[ pp. 385-387], each profile Fi is obtained 
by the minimax variational construction (see details in p4| § 3]) on the sets of the category 
p>l + l. 

Remark: convergence to periodic solutions. It is clear from two Figures above (cf. 
the last boldface profiles) that, for large I S> 1 (actually, already for / > 3 in both cases), 
the solutions of the problem Fi{y) of fl4.2p are close to a periodic structure. Namely, 
denote by F^:{y) a T^,-periodic solution of the ODE fl4.2p in M normalized so that 

sup \F^{y) \ = 1. 

Then, by scaling invariance of the equation, 

a~'^F^ (I) is also a solution in M for any a > 0. 

Therefore, for large / ^ 1, the following holds: 

(4.3) FKi/)^ar'"F,(^), where a, = f , 

and the convergence as / — > oo is uniform in y. Note that the periodic solution F^, does 
not satisfy the Dirichlet boundary conditions at y = ±R, and this creates some boundary 
layers. It is easy to see that these are of order o(l), i.e., negligible as / — )■ oo. 
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-F'®'+F'^=0, Dirichlet: first five solutions 




Figure 12. The first five solutions Fi{y) of (fO|) for m = 4. 

4.3. Homotopic connections to the cubic equation. We now introduce the basic 
one-parametric family of Dirichlet problems in {—R,R) with the operators 



n 



0, 



where e G [0, 1]. For £ = 0, we have the original problem (11. 5p . while e = 1 gives the 
above simpler problem (14. 2 p with all the solutions ordered by Sturm's index. Notice that, 
for all e G (0, 1], operators (14. 4p contain analytic nonlinearities, and the dependence on e 
is also analytic. The only problem of concern is the singular limit e — 0. 

Our further construction is naturally related to classic theory of homotopy of compact 
continuous vector fields, f2T[ § 19]. Denoting by K{y,C,) the symmetric kernel of the linear 
operator (-l)'"^^^'^) with zero Dirichlet conditions, the problems (14. 4p can be written in 
the equivalent integral form 



(4.5) ^eiF) = F-{l-e) K{F-\e^ + F^\ ■'("+'^ F) - e K F' 



0, 



where each integral Hammerstein operator is compact and continuous in (or suitable 
LP spaces for p > 2), [21, p. 83]. Therefore, the function (14. 5 p for e G [0, 1] establishes a 
deformation of the original vector field 



K{F-\f\ into <l>i{F) = F - KF'^; 



see [m p. 92]. If this deformation is non-singular (0 ^ (7($^(F))), the two vector fields 
are homotopic. For convenience, later on we consider the differential form of deformations 
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bearing in mind the necessity to return to the corresponding compact vector fields for any 
rigorous justification. 

Thus, we take an arbitrary solution Fijj) of f ll.Sp with connected symmetric (always 
achieved by shifting) compact support 

suppF = [-?/o,yo]- 

Definition 4.1. We say that a solution F{y) of (1 1 . 51) has Sturm index I = Is{F), if there 
exists its continuous non-singular deformation, also called a homotopic connection, 

(4.6) {F{y;e) on [-y,,y,]; e E [0,1]} 

consisting of critical points (solutions) of the functional for operators i \4-4^ such that 
F{y; 0) = F{y) and F{y; 1) coincides with the solution Fi{y) of (|^.^p. 



If, for a given solution F of (11. 5p . such a non-singular homotopic deformation does not 
exist, then we say that Sturm's index Is{F) cannot be attribute to such a function F in 
principle. In what follows, this nonexistence result can be associated with the fact that, 
for these solutions, the homotopic connections such as (14.51) (or others) become singular 
at some saddle-node-type (s-n) bifurcation point Ss-n < 1, at which two e-branches of 
geometrically similar solutions meet each other. 

In general, Sturm's index can be extended from e = 1 (the ordered cubic problem (14.21) ) 
to e = (the original one) along any continuous analytic branch that can have arbitrary 
even number of turning s-n points for e G (0, 1), and even beyond that. Therefore, in fact 
we ascribe the same Sturm index / to all profiles belonging to the same analytic branch 
started for e = 1 at the point Fi and ended up at e = 0. In this sense, the nonexistence 
then means that such a branch in principle is non-extensible to e = 0. 

The possibility of bifurcation (branching) points is a key difference between second and 
higher-order equations. Indeed, for m = 1, all the (or most interesting) solutions have 
the index by Sturm's Theorem on zero sets, while, for m > 2, there are many others, 
which principally cannot obey such a simple classification, associated with second-order 
problems only. 

Obviously, along any non-singular homotopic path, the critical points are deformed 
continuously, which is guaranteed by the Inverse Function Theorem {q.v. e.g. [331 P- 319]). 
Therefore, our strategy is now to use the Lusternik-Schnirel'man/fibering method for 
construction of a countable number of branches of different profiles {Fi}, which are ordered 
by the category of the sets involved. We then apply this for any e G (0, 1) in (14. 4p . Theory 
of compact integral operators [201 1211 133] then suggests existence of a countable set of 
continuous e-curves of critical points that will continuously attribute the Sturm index from 
the regular problem (14. 2 p for e = 1 io the non-Lipschitz one (II. 5p for £ = 0, provided 
that these branches are extensible and some of these are not destroyed at saddle-node (or 
others, even harder) bifurcations in between. Lusternik-Schnirel'man and fibering theory 
guarantee existence of a countable number of extensible branches. Notice that, ( see 
p. 387]), to the authors knowledge. 



"It is not known whether the Lusternik-Schnirel'man critical values are stable." 
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On the other hand, there are some definitely stable branches. Therefore, in general, we 
cannot guarantee that all the Lusternik-Schnirel'man branches are extensible to e = 1. 
In addition, the proof of the fact that the homotopic path (14.41) (or suitable others) is 
non-singular, is also a difficult open problem. 

Nevertheless, we expect that the stability or non- singularity for (14. 4p takes place for 
our particular problem, and we end up this discussion as follows: 

Conjecture |4ll. Each function Fi from the basic family {Fi, I > 0} for i\1.5\i can be 
continuously deformed by (4^) to the corresponding solutions Fi of ( |^ . . 



During the course of the inverse e-deformation from e = 1 to £ = 0, the profiles Fi{y; e) 
get a finite, depending on e G (0, 1), number of oscillations and zeros close to end points 
y = ±R, and only eventually, at e = 0, this number gets infinity, when the nonlinearity 
becomes non-Lipschitz and the solutions become compactly supported in {—R,R). 

Existence of homotopic connections for basic patterns Fi . As a typical ex- 
ample, in Figure [131 "we present such e- deformations (14. 4p of two profiles, Fo{y) (a) and 
the dipole Fi{y) (b) for n = 1 and R = 10. The e-deformation of F2{y) is presented in 
Figure ITW a) for n = 1 and R = 10. A typical corresponding e-branch of F2 is shown 
in Figure ITW b). All ^-branches of the basic family {Fi} look quite similarly. Note that 
these branches can be extended beyond e = 1, and then we observe there the absolute 
minimum of ||-F||oo at this value e = 1. 

It is an obvious observation that, by continuity of branches with respect to small changes 
of nonlinearities, if the homotopic connection as in Figure [TW b) takes place for the basic 
deformation (14. 4 p and the branch is infinitely extensible for e > 1, a. similar connection e : 
— 7- 1 can be achieved by other analytic deformations. In this sense, the type of reasonable 
homotopic deformations is not that crucial; cf. Proposition 14.21 below establishing an 
analogous non-homotopy conclusion. 

Nonexistence of ^-connections and saddle-node bifurcations . We now deal 
with other families of non-basic patterns obtained in Section [2] by an extra preliminary 
Cartesian /^-approximation. We then introduce their total generalized Sturm index, which 
should include the number of oscillations about the non-trivial equilibria ±1 defined by 
the structure of h{y); see Section [7] for an alternative approach to the generalized index 
via a spatial /^-compression of profiles. 

Then, a homotopic e-deformation of these patterns to those of the equation (14.21) with 
monotone nonlinearity is not possible in principle. We claim that, on the {e, F}-plane of 
the global bifurcation diagram, their solution branches appear in standard saddle-node 
bifurcations that occur for e < 1, i.e., these branches do not admit extensions up to the 
simpler ODEs f Oj) . 

Using the enhanced numerics with Tols=10~^ and the step Ae = 10~^, we show the 
e-deformation of two non-basic profiles given in Figure [151 

(4.7) F+^{y) and F+2,2,+2(2/), 

which have similar geometric shapes, with the equal numbers of four intersections with 
the equilibrium -|-1. This detailed e-deformation via (14. 4p for n = 1, of F+4(?/) and 
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(a) e-deformation of Fq, R = 10 
m=2, n=1 : E-deformation of F, 




(b) e-deformation of Fi, i? = 10 
Figure 13. e-deformation via (jO]) of Foiy) (a) and the 1-dipole profile Fi{y) (b). 



F2 = F+2,2,+2(y) is shown in Figure 15(a) and (b) in [13] (these are too big, 5.3 and 2.33 
MB, to be presented in arXiv). The e-deformations of these two profiles turns out to stop 
at the same saddle- node bifurcation at e = £s-n = 0.709... . 
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m=2, n=1 : e-deformation of F^^ i _2 i +2(y) 




(a) e-deformation of F2, i? = 16 
m=2, n=1: E-branch of 




(b) £-branch of F2 

Figure 14. The e-deformation of F2 via homotopy (j4.4p for n = 1 (a), and the 
e-branch of F2 (b). 



In Figure [TBI we show the corresponding e-bifurcation diagram with a saddle-node 
bifurcation at 



£s-„ = 0.709... (F+4, i?=14), 
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1.6 
1.4 



1.2 

F(y)i 

0.8 
0.6 
0.4 
0.2 


-15 -10 -5 5 10 15 

y 

Figure 15. Patterns F+4 and have similar geometric shapes and ap- 

pear simultaneously at the s-n bifurcation at (j4.8p : m = 2, N = n = 1. 

at which the branch of F+4 and the branch of F+2,2,+2 nieet each other. For convenience, 
in Figure [161 we also draw neighbouring global branches of the basic patterns F4 and 
F2 = F+2,_2,+2 (existing for all e G [0,1]), to which the corresponding branches jump 
being extended above the bifurcation values (14. 8p . 

It turns out that a neighbouring branch of the basic pattern F4 exists for e > Ss-n, while 
the neighbouring basic branch is that of F3 = F+2, 1,-2,1, +2- Being extended numerically 
for e > Es-n, the e-branches of profiles fl4.7p jump to these basic e-branches. 

Such a branching at e = Eg-n means that the two profiles fl4.7p belong to the same 
family, both having the generalized Sturm index = +4. 

It is not difficult to choose other pairs of patterns F with similar geometries, which 
have to be originated at saddle-node bifurcations for e < 1. For instance, these are F+2k 
and -F+2,2,+2,2,...,2,+2, with k single patterns ~ +Fo gluing together. 

For example. Figure [T71 shows the e- deformation of (see Figure [TST a)) 

(4-9) F+e{y) and -F+2,2,+2,2,+2(z/)- 

In Figure [TST b). we show the corresponding bifurcation diagram. Notice that the corre- 
sponding s-n bifurcation point, 

7s-n = 0.700... (F+6, R = 20)~ 

is rather close to (14. 8 p for F+4 (notice different lengths R). For F+g and -F+2, 2+2,2, +2, 2+2,2 
shown in Figure [T^ for the same i? = 20, it is different, 

£s-n = 0.52... (F+8, R = 20). 
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Figure 16. The e-branches of -F+4 and cannot be extended from e = 

via (j4.4p to e = 1 and meet at a saddle-node bifurcation at ffs-n = 0.709... . The 
neighbouring branches of -F4 and i^2j which are detected in Figure 15(a), (b) in 
|13j are also shown. 

Meanwhile, for convenience, we present the following simple conclusion showing that 
nonexistence of the homotopic path (14. 5 p actually means nonexistence of any analytic 
non-singular connections. In particular, this indicates that the geometric type of the 
branching in Figure [16] is generic. 

Proposition 4.2. Let, for a given solution F of (li.5j) . the basic deformation ( [^^.51 ) have a 
singular point at some Ss-n < 1; where two continuous branches of two patterns originated 
at e = meet each other and hence cannot be continued up to e = 1 . Then any analytic 
deformation of these patterns generating the functional path l \4-6^ ends up at a singular 
point for some e G (0, 1). 

In other words, other analytic deformations cannot move the s-n point into the set 
{e > 1} just by continuity. 

Proof. Without loss of generality, we assume that the corresponding critical values and 
the points {Fi} of the cubic problem (14.21) are non-singular (by changing R if necessary). 
Since any continuous deformation of the basic path will continuously (and analytically) 
deform the branches, the existence of a homotopic path would mean that at some instant, 
the s-n point of the branches will touch the vertical line e = 1. At this moment, we would 
create a singular value for the analytic cubic problem (14. 2p . a contradiction. □ 
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m=2, n=1 : e-deformation of F 



1.6 




y 



(b) e-deformation of i^+2,2,+2,2,+2, i? = 20 

Figure 17. e-deformation via (j4.4p for m = 2, n = 1, of F^Q{y) (a) and 
-f+2, 2,-1-2,2, +2 (y) stops at the same saddle- node bifurcation at e = es_n = 0.700... . 

Numerically, we have observed a curious phenomenon: in a left-hand neighbourhood of 
the saddle-node bifurcation at e = Ss-n, the profiles keep only essential non-monotonicity 
features and loose all intersections with zero, so become non-oscillatory near transversal 
zeros. Therefore, according to such an e-deformation to saddle-node bifurcations, the 
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0.5 0.700 1 



(b) e-bifurcation diagram 

Figure 18. m = 2, n = l: profiles -F+6(y) and F_|_2,2,+2,2,+2(y) (a); bifurcation 
diagram (b). 
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m=2, n=1 : £-deformation of F 

1.6 
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Figure 19. e-deformation via (j4.4p for m = 2, n = 1, of F^s{y) is possible 
until £s-n = 0.52.... 

number of intersections with the trivial equilibrium at e = should not be taken into 
account in the generalized Strum index. In this sense, the complicated profiles on Figure 
[TOl have the following generalized Sturm index: 

^min = {-8,+4,-10,+8-12}; 

see Section [7] for details and more mathematics. One can "split" this index to get equiv- 
alent pairs of profiles originated at some ^s-n < 1- 

4.4. Homotopic connection to linear eigenvalue problems. This is an alternative 
way to ascribe Sturm's index to basic patterns {Fi}. Modifying the nonlinearity in the 
approximation fl4.4l) . we consider the following operator family (watch the last term): 

(4.10) A,{F) = {-1)"'+^F^'^"''^ + {l-e){F-\e^ + F^\~^(^\F\'F), e G [0, 1]. 

As usual, by the actual homotopic connection we mean the corresponding vector fields 
composed of compact integral operators. Then, from (14.101) at £ = 1, we obtain the linear 
operator 

Ai(F) = (_l)"^+lF(2rn), 

On the other hand, for any e > 0, the linearized operator at is very simple, 

(4.11) a;(0) = (-1)'"+1d2'" + (1-£)/, £6(0,1]. 

Denoting by a = {— A; > 0, / > 0} the eigenvalues of the negative operator {—l)"^Dy"^, it 
follows that Ch. 8] 

El = 1 — Xi for any / > 
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are subcritical bifurcation points, where the necessary ^-branches appear. Along those 
branches that are originated at > 0, for e = we obtain our nonlinear eigenfunctions 
{Fi}, which thus inherit the Sturmian structure and the index from the eigenfunction ipi 
(it has precisely / zeros and / + 1 extrema points; see Section [6]) that governs the pattern 
for e ^ eY ■ We continue developing such an e-deformation approach to linear eigenvalue 
problems in Section |6] devoted to ODEs with analytic nonlinearities. 



5. Problem with "fast diffusion" : EXTINCTION AND BLOW-UP PHENOMENON IN 

THE DiRICHLET SETTING 

Here, using typical concavity-like techniques, we prove that finite-time extinction for 
the PDE fll.l3p is a generic property of their bounded weak solutions. Firstly, we study 
this phenomenon in a bounded domain in with homogeneous Dirichlet boundary 
conditions. Secondly, for the Cauchy problem, possible (and rather complicated) types of 
extinction patterns will be revealed in the next Section |6] by using our separate variable 
similarity patterns. 



5.1. Extinction for some nonlinear parabolic problems of higher order: main 
result. Let f2 be a bounded sufficiently smooth domain in M^. Taking the original 
equation (I1.13P and setting, as usual, 

I I — ^ 

\u\^u = V =^ u = \v\ "+1 V, where — > for n G (—1,0), 
we arrive at the following initial boundary value problem: 



+ 5 



dt 

(5.1) { v = Dv = ... = D'^-^v = on 9fi X M 

f (x, 0) = vq{x) in f2, 

where vq is an initial function from an appropriate space to be specified. Here, the only 
nonlinearity is ip{v) = \v\~'^v. We examine the problem fl5.ip in the "native" energy 
Sobolev space. 

Let us introduce the following functionals associated with the operators in (15. ip : 



*(^) ■= I / \'v{t,x)\"+^^ dx, 
Jn 

E{t) :=- [ \D"'v{t,x)\^dx- I v'^dx. 
Jn Jn 



Lemma 5.1. There holds 

$'(t) = n±lE{t). 

The proof follows from simple calculations. The main result on extinction in (15. ip is as 
follows: 
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Theorem 5.1. (Extinction) For given nontrivial initial data, denote: 

$(0)= i / \vo\^dx > 0, 
Jn 



E{0) = - I iD^'vol^dx- I vldx<0. 



Then 

I "+2 I 

(5.2) $(t) < $(0) = 1 - ^) ' " ' ^(0) ^0 as t^T-, 

(i-i) " 

with 

rn ._ n+2 ^{0) _ 2 ^{0) ^ r, 
• n ^'{O) n^;{0) ^ 

The proof is divided into several steps. 

5.2. The first energy relation. Multiplying the equation in (15.11) by v and integrating 
by parts over VL x (0,T) by taking into account the boundary conditions, we obtain 



(5.3) 



1 

n+2 



n + 2 

\v\ "+1 



dx = — 



D"'v\'^ dxdt 



V dx dt. 



5.3. The second energy relation. Multiplying the equation in (15.11) by vt and again 
integrating by parts over Q x (0,T) by using the boundary conditions, we obtain 



(5.4) 



%()'{v)v^ dx dt 



iD'^v^dx 



V dx 



5.4. Connection (a main point). Denote 

U{t):=- [ iD'^vl^dx- I |t;pdx, V{t) := ^ I |t;|STfdx. 

Jn " 

Then the identity (15. 3p reads 



(5.5) 



Vit) 



T 



dV 

U{t)dt, or equivalent ly, -— = U(t). 



dt 



Analogously, it follows from (15. 4p that 



(5.6) 



dU _ 2 
dt ~ n + 1 



\v\ "+i (t>t)^ dx. 
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5.5. The main (crucial) relation. We deduce from (15. 5p and (15.61) that 
(5.7) d^=d^ = ;^/l^l 



n 

Next, by the definition of V , we have 

[n + I) Jq^ n + i Jq 

Thus, equation (15. 7p takes the form 

(5.8) / \v\~^vvttdx = / \v\~^{vt)'^ dx. 

Jn n + 1 Jq 

5.6. Replacement. We introduce a new function z = z{x,t) by the formula 

n 

V = \z\"+'2Z. 

Then equation (15. 8 p reads 

(5.9) / zztt dx = Cn / {ztY dx, with c„ 



n Jn 



,2 , 3n + 2 



5.7. Fourier analysis of equation (15. 9p . Let {ek)ke'N be a complete orthonormal system 
in L^(f2). Then we have 

z{x,t) = T.kLi^k{t)ek{x), 

ztt{x,t) = J2T=i4itMx). 

Let 



If 1 °^ 

m = 7^ / z'{x,t)dx = -J2^k{tY, 

^ k=i 



so that 

(5.10) $'(t) = 

^ ^ ^ ^ 2 dt 



« « CX) 

/ dx = / zzt dx = 2_\ Zkz'k, 
Jn Jn f.^^ 



and hence, by Holder's inequality, 

(oo \ / oo \ oo 
k=l ) \k=\ J k=l 

In Fourier coefficients, equation (15.91) takes the form 

oo oo 

(5.11) ^z4' = c„J]z^ 

k=l k=l 
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Zk ■ 



On the other hand, in terms of the function $, we have 

oo oo 
fc=l k=l 

Consequently, equation (15.111) takes the form 

oo oo 



k=l k=l 



Therefore, 



k = l 

12 



I.e. 



(5.12) > where 

fc„ _i±^^ 2(n + l) (forne (-2,0)). 

5.8. Extinction: the proof. For the analysis of an ordinary differential inequality ap- 
peared, we can use various approaches. For the case where $' < 0, we apply the standard 
approach. Namely, we divide inequality fl5.12p by Then we obtain 

$" $' 

From here, it then follows that 

$'(t) < Ci$^"(t), 

$'(0) 

with Ci = \ < for E(0) < 0. 
$'^"(0) 

This inequality implies the result of the main theorem. For the proof, it suffices to 
replace our new variables with the original function v by the formula 



so that 



$(t) = - / z'^{x,t)dx = - / dx 
2 in Jn 



N f . 1^ + 2/",,, , 

$ (t) = — = / zzt dx = - — — - / tp{v)vt dx 
dt J<^ 2 n+1 J<^ 

1 n + 2 f il){v 



2n+l^ ij'{v 

ln + 2 
' 2 n+1 

n + 2 
' 2 

n + 2 
2 
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ip'{v)vt dx 
n + l)vilj{v)t dx 
vip{v)t dx 



E{v), 



where 

E{v) = - I dx- [ dx. 

Jn Jn 
This completes the proof of Theorem 15. 1[ 

5.9. A diversion to blow-up for n > 0. By using the same computations as above we 
can prove the following blow-up result for the original equation (II. ip . 

Theorem 5.2. (Blow-up) Suppose that n > in fl5.jp . Let 



If n + 2 

$(0) = - / |fo|"+i dx > 0, 
2 Jq 

E{0) = - [ \D"'vo\^dx+ [ v^dx>0. 
Jn Jn 



Then 

m > m- — ^—^ 

{l~t/T)~ 

with 

n + 2$(0) 2<l>(0) 



n $'(0) nE{{])' 



5.10. Blow-up: the proof. In this case, following the computations as in (15. lip , we 
obtain, 



Therefore, 



k=l k=l 
k=l 

i.e., 

(5.13) > 
where 

1 + Cn n + 1 
K := = 2^— > 1 for n > . 

2 n+2 ^ ' 

For the case where $' > 0, we apply the standard approach. 

Namely, we divide inequality (15.130 by Then we obtain 

$" $' 

From here, it follows that 

$'(t) > Ci$(t)^" 

$^0) 
$(0)* 
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with Ci = —\t > for E(0) > 0. 



This inequality implies the result of the main theorem. 

Indeed, it is enough to replace our new variables with the original function v by the 
formula 

<l>(t) = - / z(t,xf dx = - \v\^ dx 



2 



n 



and 



N d$ f , 1 n + 2 f 

$' t = — = / zztdx = / 

dt Jn 2 n + lU 



ip{y)vt dx 



ip [v)Vt dx 



2 n + lJnij'{v) 
1 n + 2 



2 n + lJn 
n + 2 
2 

n + 2 



n + l)vilj{v)t dx 
vtp{y)t dx 



n 



2 



E(v] 



where 

E(v) = - f ( \D"'v\'^ + v^] dx. 



6. Problem "fast diffusion": EXISTENCE, MULTIPLICITY, ETC. 

6.1. Oscillatory ODEs with analytic nonlinearities from fast diffusion. Here we 
consider another ODE model (11. 7p . and without loss of generality, we mainly restrict to 
m = 2. We show that (ll.7p provides us with similar countable families of various patterns. 
Moreover, we claim that the solution set of (II. 7p is equivalent to that obtained earlier for 
non-Lipschitz nonlinearities. Indeed, solutions F{y) of (ll.7p are not compactly supported 
and exhibit oscillatory exponential decay at infinity governed by the linearized operator 

(6.1) F(^) = -F + ... =^ F{y) ^ e~^cos{^y + c), y ^ +oo. 

6.2. Patterns. Figures |20] and |2T] show a few typical patterns, which we are already 
familiar with. It is important to notice that, in Figure [201 by watching the behaviour of 
small negative solutions close to the origin y = 0, there are two different patterns that can 
be classified as F+2,2,+2, and the second one is denoted by F^^ 2 +2- This shows again that 
the number of intersections with equilibria ±1 and are not enough for a complete pattern 
description (in fact, this underlines that a homotopy approach using the hodograph plane 
is not applicable to equations such as (11.160 ). It is seen there that -^+2,2, +2(2/) exhibit more 
non-monotone structure for y G (—4, 4) than F+2,2,+2(|/), so that the derivative F^2 2 +2(2/) 
has there 3 zeros therein, while -F+2,2,+2(2/) has just one at y = 0. Thus, these two patterns 
can be distinguished by the number of zeros of their derivatives, but, by the same reasons, 
we do not think that counting internal zeros of the pairs {F{y), F'{y)} can help to create 
any rigorous Sturmian-like classification of such patterns. 
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y 



Figure 20. Various solutions of the ODE (jl.7p : positively dominant patterns. 

In Figure [221 we present typical complicated patterns for the problem (11.7p . which 
remind similar "multi-hump" structures obtained above in Figure [THl 

Note that, in view of the fast exponential decay (16. ip . it is difficult to observe by 
standard numerical methods that, unlike the previous problem, the profiles in Figure [22] 
are not compactly supported. Note that, for n > 0, we succeeded in the logarithmic scale, 
by taking small regularization parameter and tolerances ~ 10~^^, to reveal the difference 
between the linearized zero-behaviour (as in (16. ip ) and the nonlinear one at the interface 
as y ^ Uq of the type {f{s) is the oscillatory component; see Figure [6] and [HI [121 § 4]) 

(6.2) F{y) = {y- yo)"'(f{s), where s = ln{y - y^) and 7 = Mn±li > 2m. 

For n < 0, there are no "nonlinear zeros", but there exists an infinite nnm^er of linearized 
ones given by (16. ip . first dozens of which can be easily observed numerically. 

6.3. Application of Lusternik— Schnirel'man and fibering theory. Obviously, the 
ODE (11.70 (or the elliptic problem in (I1.16P ) possesses the variational setting in {—R,R) 
(or Bn), to which the same fibering version of Lusternik-Schnirel'man theory applies, as 
in fn\ [T21 § 3]. This gives a countable family of basic patterns {Fi} for both the ODE 
and the elliptic PDF. Introducing the preliminary /i-approximation of patterns makes it 
possible to reconstruct other families of more complicated geometric structure. Since Fi 
are not compactly supported, we always assume fixing R = R{1) 1 for / 1. The 
solutions in M (M^) are then obtained by passing to the limit / — > 00. We do not stress the 
attention to such a compactness procedure that assumes deriving some uniform bounds 
independent of /. On the other hand, in view of the known exponential decay of all the 
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m=2, N=1 : F<'*'=-F+F^ 




Figure 21. Various solutions of the ODE ()1.7|) : essentially changing sign patterns. 

solutions, the variational statement in the whole space is also an option; see [26] . We 
now present brief comments. 
The functional is 



(6.3) 



The spherical fibering [TTl [121 § 3] with v belonging to the unit sphere, 
(6.4) F = r{v)v, where ve'Ho = {Ho{v) = j \E)'^v\^ + / " 

leads to the functional 
(6.5) 

This attains the absolute maximum at 
(6.6) Hl = r-r^ I = 



H{r,v) = \ r'^ — \ r'^ I v'^. 



ro(v 



at which H{ro{v),v) 



(6.7) 



= — . This defines the positive homogeneous convex functional 



-^('^) ~ [2H(ro{v),v)] 
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F(y) 



-0.5 - 



p(4)_p3_p^ Dirichlet: a complex pattern 





-60 -40 -20 



20 40 60 



(a)F, 



+ 10,l:-4,+6,4,+8 




100 



(b) +4. 1,-12,1, + 12, 1,-16 

Figure 22. Examples of complicated patterns of the ODE (jl.7p . 



Here Lusternik-Schnirel'man theory applies in its classic form |2T1 p. 387] giving a count- 
able set {cfc} of critical values and points denoted by {vk} (see full details in [HI [121 
§3]): 



(6.8) 



Cfc = supj-g_v(fc inf^eJ- j v'^- 
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Now the sequence of critical values is decreasing, 



(6.9) ci > C2 > ... > Cfc > Cfc+i > ... . 
For each Vk (a solution F^), the critical values are given by 

(6.10) v = CF eUo =^ cf = H{v) = 



For all the patterns shown in Figures [20] and |5T1 these values (Lusternik-Schnirel'man 
critical or not) are presented in Table 1. As usual, this table makes it possible to detect 
the Lusternik-Schnirel'man critical points that deliver critical values (16. 8p for each genus. 
We then again claim that the basic patterns {FJ} deliver all the Lusternik-Schnirel'man 
critical values (16. 8 p with k = 1 W. 

Recall the Formal Rule of Composition (FRC) from pTl [T2l § 5] : performing maximiza- 
tion of H{v) of any {k — l)-dimensional manifold J-" G Mk, 
(6.11) 

the Lusternik-Schnirel'man point Fk-i{y) is obtained by minimizing all internal tails and 

i.e., making the minimal number of internal transversal zeros between single structures. 
Then (I6.1ip also applies, since by the same reason diminishing a small tail between two 
Fo-structures will increase the corresponding value cp in (I6.10p . 



Table 1. Critical values of H{v) = 




F 




Cf 




Fo 




0.2033.. 


= Ci 


Fi 




0.1080.. 


= C2 


F+2,2,- 


h2 


0.1019.. 




F+2,2,- 


h2 


0.1017.. 




F-2,1- 


h2 


0.1014.. 




F+, 




0.0961.. 




F2 




0.0736.. 


= C3 


F+2,2,- 


h2,2,+2 


0.0680.. 




F+2,1,- 


-2,1, +2 


0.0675.. 




F+e 




0.0629.. 
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6.4. Homotopic connections with Sturm-ordered hnear eigenvalue problems. 

An extra advantage of the problem fll.7l) is that the homotopic connections of basic 
patterns can be revealed easier. Namely, as above, we consider the ODE (11. 7p in {—R, R) 
with sufficiently large R > to see the first / patterns which are still exponentially small 
for y ^ ±R. Consider the following homotopic path with the operators: 

(6.12) A,{F) = F^^^+eF-F^ = 0. 
Consider the corresponding linearized operator: 

(6.13) A'M = D'y + eL 

Let cr{Dy) = {A/>0,/>0}be the discrete spectrum of simple eigenvalues of Dy > 
in L'^{{—R,R)) with the Dirichlet boundary conditions. The orthonormal eigenf unctions 
{ipi} satisfy Sturm's zero property; we again refer to [7j for most general results. 
By classic bifurcation theory [211 P- 391], for such variational problems, 

5; = -A, < for / = 0,1,2,... 

are bifurcation points, so there exists a countable number of branches emanating from 
these points (but we take into account the first ones). In order to identify the type of 
bifurcations, in a standard manner, setting e = ei + s for |s| small and 

F = Cipi + w, where W-Lipi, 

substituting into fl6.12p and multiplying by ipi yields 



C\s) H^f + ...>0 ^ s> and C{s) = ± 



\ 



+ 



Hence, at e = ei, there appear two branches with the equations 

Fi{y) = ± l^f^ i,i{y) + ... for e > Si. 




Therefore, all the bifurcations are pitchfork and are supercritical, i.e., two symmetric 
branches are initiated aX e = see Figure I^ST a). where we construct numerically the 
first positive e-branch of Fq of the e-bifurcation diagram and show that this branch is 
extensible up to the necessary value e = 1, and even up to e ~ 10^ (b), and further. 

It follows that all the branches are originated at < 0, so being continued up to e = 1, 
give the original equation (11. 7p . The questions of global continuation of branches are 
classic in nonlinear variational theory; see P, § 6.7C]. The global behaviour of bifurcation 
branches for 2mth-order ODEs with analysis of possible types of end-points is addressed in 
[2]. These results hardly apply to the equation (14. 4p with non-coercive operators admitting 
solutions of changing sign near boundary points. 

The existence of a turning point of the given branch in this real self-adjoint case, i.e., 
of a saddle-node bifurcation, assumes that there exists an eigenvalue (say, simple) 

Oea(A;(F)), i.e., 3 0o: 0^'^ + £0o - SF^o = 0, 
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where (po is an eigenfunction of A^(F) satisfying the Dirichlet conditions (po = 0o = at 
y = ±R. For a moment, we digress from our difficult ODEs and consider simpler models 
with known bifurcation diagrams. 

Remark 1: on turning saddle-node bifurcation points of positive solutions. 

Such turning points do exist for equations with other nonlinearities and another depen- 
dence on e, e.g., 

(6.14) = ee^ or (-l)™^^^"^) = e(l + (F > 0), 

with Dirichlet boundary conditions at ±R. Existence of two different branches of solutions 
(they are positive by [7]) for all small e > is established by the fibering method as in 

[mini §3]. 

On the other hand, for £ ^ 1, positive solutions of f l6.14p are obviously not possible. 
This is easily seen by multiplying the first equation in ( I6.14p by the first eigenfunction 
ipo > with the eigenvalue Aq > of the positive operator (— l)™/)^™. This gives for the 



first Fourier coefficient Cq = j Fipo the following inequality: 



(6.15) \oCo = e J e'\p,^e\mij e^j^^>e\mieJ ^H^olli = e||^o||i eW^ft. 

At the last stage we have used Jensen's inequality for the convex function e^. One can 
see that the resulting inequality in f l6.15p . 



does not have a solution Co > for all £ ^ 1, meaning nonexistence. 
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14.91 



Figure 24. The bifurcation diagram for (j6.16p . 

By standard theory of compact integral Uryson-Hammerstein operators [201 EI] and 
branching theory |33], the solutions of f l6.14p detected for small e > comprise two 
continuous branches, which must end up at a saddle-node bifurcation point (no other 
bifurcations are possible) at some 

e = es-n > 0. 

In Figure [2H the global bifurcation diagram is presented for the quadratic equation, 
(6.16) F(^)=e(l + F2) on (-1,1), where 



14.91... 



The upper branch blows-up to +oo as e — )■ 0"^, while the lower one vanishes. See [19] for 
a different approach to bifurcation analysis for such fourth-order ODEs. 

We return to the equation f l6.12p . Nonexistence of such saddle-node bifurcation points 
for the ODEs fl6.12p is a difficult open problem. Numerically, we have got a strong evidence 
that the first branches are strictly monotone without turning points; see Figure IMT b). 
where the first branch is extended up to e = 100. 

On the other hand, it is easy to check by the fibering approach that f l6.12l) has a 
countable set of solutions for arbitrarily large e > 0, which are also identified as continuous 
curves by nonlinear compact operator theory. Hence, there are infinitely many branches 
of basic patterns that are unboundently extensible in e. We then conclude that at e = 1 
all these branches are available. Hence, at e = 1, the corresponding solutions Fi of (11.71) 
inherit their Sturm index from eigenfunctions ipi that occur at e = e^. 
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Of course, we can use the alternative homotopy approach. Figure \23( a) shows that the 
e-branches are well-defined at e = 0, where we obtain the simpler equation 

that, as in (14. 2p . admits Sturm's classification of all the patterns. 

Remark 2: explicit global monotone ^-bifurcation branches for a non-local 
nonlinearity. Computations of branches are straightforward for the following non-local 
equation with the analytic cubic nonlinearity (cf. (I6.12p ): 

(6.17) F^^^ = -eF + {J F^)F in {-R,R). 

Then the solutions 

Fz(y) = ±v/A7Ti V^z(y) for / = 0,1,2,... 
are originated at ei = — A^^, and the branches have the geometric form as in Figure [251 

7. Problem "Sturm Index": USING _R-COMPRESSION 
Consider again for a while the simplest non-Lipschitz ODE for m = 2: 

(7.1) F(^) =F -\f\ in M. 

As we have seen, the oscillatory character of the solutions close to interfaces, finite or 
infinite, causes some difficulties in determining the actual (generalized) Sturm index of 
the patterns defined as the number of certain oscillations, which are "dominant" zeros 
(or extrema) points about 0. Oscillations about non-trivial equilibria ±1 are clear and 
do not exhibit such difficulties, though sometimes it is rather hard to distinguish these 
classes of non-monotonicities. The main problem of concern is that how to distinguish 
the "dominant" non-monotonicities and "small" non-dominant ones related to asymptotic 
oscillations close to interfaces, which are almost the same for all the patterns; cf. [HI [T21 
§ 4] for the non-Lipschitz nonlinearity and (16. ip for the analytic one (II. 7p . 

As we know, in the CP for the PDEs involved, it is enough to pose the Dirichlet problem 
for (ll.Sp on sufficiently large intervals {—R, R), so one does not need to take infinite length 
R = +00. Let is concentrate on the case m = 2, 

(7.2) FW = F-|F| in {-R,R), 

with Dirichlet conditions at the end-points y = ±R. Given a compactly supported pattern 

(7.3) F{y) = Ftj{y), with multiindex a = {cxi, (T2, a^, cxi 

we perform its R- compression, i.e., start to decrease R observing a continuous deformation 
of the corresponding profile F{y; R) until the minimally possible value 

R = Rmin > 2i?, > {R^ is defined in [HI [121 § 3] ) . 

Let (Tmin be the multiindex of F{y; R^ij^), if this profile is bounded. If not, we mean 
the index F{y; R^^^) calculated for profiles with R ^ -Rmin- usual, (Tmin refiects the 
sequence of intersection numbers with equilibria ±1, without taking into account those 
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with 0, which are, actually, nonexistent; see below. There is a direct relation to Definition 
l4.1l to be explained later on: 



Definition 7.1. Given a compactly supported solution ( |y.g| ) of ( | 7. , by its generalized 



Sturm index we mean the multiindex o^am obtained by the R- compression. 
Two cases of i?-compression are distinguished: 

(i) As -R — )■ -Rmin) profile F{y] R) gets unbounded and achieves a clear geometric 
form with / extrema. Then we say that this / is precisely the Sturm index / = Is of the 
functions F on this i?-branch. As usual, we claim that such an index / can be attributed 
to the basic family {Fi} only. 

(ii) There exists a finite limit F{y; Rmin)- Then the only generalized Sturm index can be 
defined as a characterization of its geometric structure with fixed numbers of intersections 
with equilibria ±1. 

In Figure [251 we present some numerical results of the /^-compression of the profile 
Fq (a) (interfaces become non-oscillatory already for i? ~ 5, cTmm = {+2} as expected, 
i.e., with no essential "transversal" zeros and a single maximum) and F+2,2,+2 (b) (non- 
oscillatory interfaces for R ~ 10, (Train = {+4}, no essential zeros and three extrema). In 
(b), by dotted fines, we denote some other profiles that also appear for such R ^ -Rmin- 
These belong to other branches of solutions of fl7.2p . 

In Figure [221 we show the /^-compression of a different profile, F+4, with a similar 
structure. Note that, in Figures [26] and [25r b). the profiles for i? = 8 coincide. This 
again confirms (cf. Figure [16] for an e-deformation of those) that these profiles belong 
to two branches originated at a supercritical saddle-node /^-bifurcation at -Rmin ~ 7.9. 
Observe in Figure [26] other dotted profiles of a similar geometric structure that exist close 
to -R^in indeed correspond to the third basic pattern F2 = F+2 ,-2,+2- The intriguing 
global /^-bifurcation diagram with saddle-node bifurcations turns out to be similar to that 
observed in Section [H by the e-homotopy approach; see more explanations below. 

We do not need to develop more consistent theory of /?- compression in view of the 
following simple comment: 

£-HOMOTOPY AND /^-COMPRESSION CAN BE EQUIVALENT . For simplicity, consider the 
analytic ODE problem fl6.12p on (—/?,/?), where we perform the scaling 

F{y) = y/^V{z), y = eh {e > 0). 
Then e is scaled out from the ODE and enters the interval, where the problem is posed: 
(7.4) V^^'^ = -V + V^ on {-eiR,e-^R). 

Therefore, in this particular case, the e-deformation of the equation f l6.12p is equivalent 
to the /^-compression for (17.40 as e > decreases. So these two approaches to Sturm's 
index of highly oscillatory structures are essentially equivalent and hence lead to the same 
results. In particular, this explains the striking phenomenon (observed in Section [1]) that, 
at saddle-node bifurcations, the profiles lose all their oscillations about zero. 
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(a) Fq 



m=2, n=1 : R-compression of F^^ 2 +2(y)' '^min ~ 




(b) 



Figure 25. The ii-compression of solutions of the ODE fTTI) for n = 1; Fq (a) 
and -F+2,2,+2 (b). 
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m=2, n=1 : R-compression of F (y), R . -7.9 



F(y) 















R=15H 


: \ R=8 







-10 





y 



10 



Figure 26. The i?-compression of the solution F+4 of the ODE ([7TT|) for n = 1. 

8. On essentially quasilinear extensions: gradient diffusivity 

8.1. Variational problems with p-Laplacian operators. For simplicity, we consider 
PDEs (11.201) in ID, where these have simpler forms. 



(I) 
(11) 
(HI) 



"xxi 'Uxx)xx + I^^Tm (parabolic), 

'^Uxx)xx + H'^u (hyperbohc), 
"xxi ^u^x)xxx + (|MrM)x (NDE). 



Ut = -(IMi 
Utt = -{\Ux 
Ut = -{\Ua 

For the reaction-diffusion PDF (I), the blow-up solutions are the same, (II. 2p . with 

(8.2) ^/ = -(irr/T + i/r/- 

Using the scaling yields the basic quasilinear ODF: 

(8.3) / 



n 



-( 


F" 




F 


''F -F = in R. 



For the hyperbolic PDF (II), we construct the blow-up patterns (II. 9p . where the same 
scaling yields (18. 3p . Finally, for the NDF (III), the TW compacton 

Uc{x,t) = f{x — Xt), with A = — 1, 



directly leads to the ODF in fl83|) . 

In all the three cases, for the A^-dimensional PDFs (ll.20p . we arrive at the elliptic PDF, 



-A(|AF 


"AF) + 


F 


"F - F = in R, 
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admitting variational formulation in L^(R^) with the potential (cf. (12. ip ) 
(8.5) E{F) = — ^ / |AFr+2 + ^L I \F\^+^ ^ 



""''^ JiRJv n + 2 J^N 2 

Then the same Lusternik-Schnirel'man/fibering approaches can be applied to get a count- 
able basic family of compactly supported solutions and next discover other countable sets 
of blow-up patterns, etc. We claim that the most of principal results obtained above 
for the semz/mear elliptic problems can be extended to the quasilinear one (18. 4p . though 
some local (e.g., oscillatory structures of solutions near finite interfaces) or global aspects 
of the behaviour of patterns are more complicated. Some local oscillatory properties of 
solutions of changing sign are discussed in [151 pp. 246-249] and [5]. All these problems 
admit more general 2mth-order extensions along the same lines as usual. 

8.2. Related non-variational problems: branching "from" potential results. Fi- 
nally, let us mention that there exists a variety of slightly changed PDEs (I)— (II), which, 
on exact blow-up solutions, reduce to elliptic or ODE problems that are principally non- 
variational. A systematic study of such problems with non-coercive, non-monotone, and 
non-potential operators touches another wide area of open mathematical problems. 
For instance, consider a reaction-diffusion equation with the nonlinearity from (ll.2ip . 

(8.6) ut = -A(|A?i|"AM) + \u\^u, where p > n > 0. 

The similarity blow-up solutions are not separable as in (II. 2p and have the form 

(8.7) ns(x,t) = (T-t)-i/(y), y = x/{T-t)^, where = > 0. 
Substituting (18. 7p into ( 18. 6p yields the following elliptic equation: 



-A(|A/rA/)-/3yV/-i/+|/|^'/ = in 



For p = n we have /3 = 0, so that (18. 8p reduces to the above variational equation (18.40 . 

We claim that, for p n, the differential operator in (18. 8p is not variational (even 
for n = 0; see |T71 § 7]), so all above techniques fail. Nevertheless, we also claim that 
the present variational analysis can and does play a role for such problems. Namely, the 
variational problem (18.41) for p = n has the following meaning: 



.9) 



potential eq. ( 18. 4 p initiates branching at p = n of solutions of non-potential 
eq. ( 18. 8p from the Lusternik-Schnirel'man/fibering patterns {Fi}. 



Such ideas well-correspond to classic branching theory; see Vainberg-Trenogin [33]. There- 
fore, we expect that still, even not being variational, 

(8.10) (18.81) admits arbitrarily large number of solutions for p ^ . 

Here we should assume that p is sufficiently close to n"*", since, as usual, global continuation 
of local bifurcation branches is a difficult open problem. Note that, for p > n, solutions 
of (18. 8p are not compactly supported in general (for p G (0,n] these are). Applications of 
such a branching approach to (18. 8p with = 1 and p > n and p < n are given in [9] . 
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Branching phenomena for such nonhnear and degenerate operators as in (18. 8p are not 
standard and demand a lot of work; see e.g., [I], where thin film operators have been 
dealt with. On the other hand, the p-branching analysis in the semilinear case n = for 
equations such as (18.81) uses spectral theory of non self-adjoint operators and is easier; see 
examples in [5l [TOl [17] . 

In other words, precisely (18.91) is the actual role that variational problems can play for 
describing finite or countable sets of solutions of principally non- variational ones, so that 
the results such as (18.101) can be inherited from a suitable potential asymptotic setting. 
Finding such a variational problem by introducing a parameter as a good approximation 
of the given non-potential one can be rather tricky, though for equations like (18. 8p . this 
looks natural and straightforward. 
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